Nonlinear level crossing models 
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We examine the effect of nonlinearity at a level crossing on the probability for nonadiabatic tran- 
sitions P. By using the Dykhne-Davis-Pechukas formula, we derive simple analytic estimates for 
P for two types of nonlinear crossings. In the first type, the nonlinearity in the detuning appears 
as a perturbative correction to the dominant linear time dependence. Then appreciable deviations 
from the Landau-Zener probability Plz are found to appear for large couplings only, when P is 
very small; this explains why the Landau-Zener model is often seen to provide more accurate re- 
sults than expected. In the second type of nonlinearity, called essential nonlinearity, the detuning 
is proportional to an odd power of time. Then the nonadiabatic probability P is qualitatively and 
quantitatively different from Plz because on the one hand, it vanishes in an oscillatory manner as 
the coupling increases, and on the other, it is much larger than Plz- We suggest an experimental 
situation when this deviation can be observed. 



PACS numbers: 03.65.Ge, 32.80.Bx, 34.70. +e, 42.50.Vk 



I. INTRODUCTION 

The Landau-Zener (LZ) model Jl|,|| is widely used in 
quantum physics to describe level crossing and avoided 
crossing transitions. It provides the probability of tran- 
sition between two quantum states coupled by an ex- 
ternal field of constant amplitude and time-dependent 
frequency which passes through resonance with the tran- 
sition frequency. This level crossing, seen in the diabatic 
basis (the basis of the two bare states), appears as an 
avoided crossing in the adiabatic basis (the basis com- 
prising the two eigenstates of the Hamiltonian). Cases 
of level crossings and avoided crossings can be met in 
many areas in physics, such as laser-atom interactions 
||, magnetic resonance pL slow || and cold Q atomic 
collisions, molecular physics [[7]||, optical atoms P,|lO|], 
^olid-state physics []L3| -flq1 , ultrasmall 
[L9[, nuclear physics |2C| ], and parti- 
The LZ model is a basic tool for 
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tunnel junctions 
cle physics 

describing and understanding such phenomena. 

Given that the LZ model presumes very crude time de- 
pendences for the coupling (constant) and the detuning 
(linear), it is somewhat surprising that it has been often 
found to provide rather accurate results when applied to 
specific cases involving time-dependent couplings (e.g., 
pulse-shaped) and nonlinear detunings. To the best of 
our knowledge, no satisfactory explanation of this fact 
has been given so far. In the present paper, we try to 
answer this question by considering models where the 
coupling is still constant but the detuning is nonlinear. 
We investigate two main classes of nonlinear detunings. 
In the first class, the nonlinearity appears as a correction 
to a dominant linear time dependence near the cross- 
ing; we call this perturbative nonlinearity. In the second 
class, the detuning is proportional to an odd power of 
time, A oc t , with N = 3, 5, 7, . . .; hence, the detuning 
cannot be linearized in any vicinity of the crossing and 
we call this nonlinearity essential. We have found sig- 



nificant qualitative differences between these two cases. 
We study the nonlinearity effects both numerically and 
analytically and we are interested in the near-adiabatic 
regime which is the one of primary interest as far as level 
crossing models are concerned. The analytic approach is 
based on the Dykhne-Davis-Pechukas formula |25|,^6| and 
a generalization of it that accounts for multiple transition 
points |26|-|0|. 

The paper is organized as follows. The basic equations 
and definitions are given in Sec. Q. The perturbatively 
nonlinear models are discussed in Sec. [II and the essen- 
tially nonlinear models in Sec. The conclusions are 
summarized in Sec. [v[ The numerical method used for 
highly accurate integration of the Schrdinger equation is 
nontrivial and it is presented in the Appendix. 



II. BACKGROUND 



A. Basic equations and definitions 



We wish to solve the Schrodinger equation (h = 1), 



ij t <t) = H(i)c(i), 



with the Hamiltonian 



H(*) = 



-A(t) Q(t) 

n(t) A(t) 



(i) 



(2) 



and c(t) — [ci(t), C2(t)] , where c\(t) and C2(t) are the 
probability amplitudes of states ipi and ^i- The detuning 
A(t) and the coupling il(t) are both assumed real. We 
are interested in the case when the coupling is constant 
and the detuning is an odd function of time, 

Q(t) = const, A(-i) = -A(t). (3) 

Such is the case for the Landau-Zener (LZ) model, 



1 



Sl(t) = const, 



(t) = ph. 



(4) 



Since the transition probability is invariant upon the sign 
inversions A — > — A and VI — > — fi, we assume for sim- 
plicity and without loss of generality that both Vt and 
the slope j3 2 of the detuning at the crossing (t — 0) are 
positive and that j3 itself is positive. 

The level crossing, seen in the diabatic basis (|l|), trans- 
lates into an avoided crossing in the adiabatic basis be- 
cause the eigenvalues ±£ (t) of the adiabatic states (the 
quasienergies) do not cross but only come close to each 
other near t — 0. Here 



£(t) = v^ 2 + A 2 (t). 



(5) 



If the system is initially in state [c\{— oo) = 1, 
ci{— oo) = 0], the transition probability in the diabatic 
basis at +oo is given by |c2(+oo)| 2 . The transition prob- 
ability in the adiabatic basis (the probability for nonadi- 
abatic transitions) is equal to the probability of no tran- 
sition in the diabatic basis, P = |ci(+oo)| , and its de- 
termination will be our main concern. For the LZ model, 



P = e 



(6) 



where a = f2//3 is the normalized coupling. 

In Sees. [II and |y| we consider two types of nonlinear 
detunings and we investigate how the nonlinearity affects 
the LZ probability (|^). In the first type, the nonlinearity 
appears as a perturbative correction to a dominant lin- 
ear time dep enden ce near the crossing; for the superlinear 
model (Sec. Ill A ) this cor rectio n is positive whereas for 
the sublinear model (Sec. IIIB) the correction is nega- 
tive. In the second type, the detuning is proportional to 
t , with N = 3, 5, 7, . . ., and the detuning is essentially 
nonlinear. Like in the LZ model, the detuning diverges as 
t — ► ±oo in all considered cases. The studied nonlinear 
detunings are shown in Fig. [l]. 



B. Dykhne-Davis-Pechukas formula 

The Dykhne-Davis-Pechukas (DDP) formula [||J26| 
provides the asymptotically exact probability for nona- 
diabatic transitions in the adiabatic limit. It reads 



-2lmV(t c ) 



where 



V(t c ) = 2 / £{t)dt. 
Jo 



(7) 



(8) 



The point t c is called the transition point and it is defined 
as the (complex) zero of the quasienergy, 



£{t c ) = 0, 



(9) 



which lies in the upper complex t- plane (i.e., Im£ c > 0) 
and if there are more than one such zero points, t c is the 
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FIG. 1. The level crossing cases studied in this paper. 
Upper figure, the detunings A(t) in the diabatic basis. Lower 
figure, the quasienergies ±£ (t) in the adiabatic basis. 

one closest to the real axis. Equation (Q) gives the cor- 
rect asymptotic probability for nonadiabatic transitions 
provided (i) the quasienergy £(t) does not vanish for real 
t, including at ±oo; (ii) £(t) is analytic and single- valued 
at least throughout a region of the complex i-plane that 
includes the region from the real axis to the transition 
point t c ; (hi) the transition point t c is well separated 
from the other quasienergy zero points (if any) and from 
possible singularities. Amazingly, for the LZ model, the 
DDP formula (@) gives the exact probability (||) not only 
in the adiabatic limit (a 3> 1) but also for any a. 

In the case of more than one zero points in the upper t- 
plane, Davis and Pechukas [[26] have suggested, following 
George and Lin [§7|], that Eq. (@) can be generalized to 
include the contributions from all these N zero points 
t k c in a coherent sum. This suggestion has been later 
verified by Suominen, Garraway and Stenholm [^8] [H). 
The generalized DDP formula has the form 



2 



p 



N 



k=l 



where 



T k =4i lim (i-i?Wi) 

t—>ti 



(10) 



(11) 



and d(t) = [fl(t)A(i) - fi(*)A(*)]/2£ 2 (t) is the nonadia- 
batic coupling, where an overdot means a time derivative. 
In principle, Eq. ( p^| ) should be used when there are more 
than one zero points which are closest to the real axis and 
have equal imaginary parts and Eq. ( [To| ) should include 
only the contributions from these zeroes. The contri- 
butions from the farther zeroes are exponentially small 
compared to the dominant ones and may therefore be ne- 
glected. Retaining the contributions from all transition 
points, however, may be beneficial and it has been shown 
elsewhere [|29|j3l|] that for the Demkov-Kunike models f3"3] 



the full summation in Eq. (10) leads to the exact result. 

The motivation for developing good approximations, 
such as Eq. ( |l0| ) for the transition probability in the adi- 
abatic limit, does not arise only from the general wish 
to have analytic expressions. In fact, the numerical in- 
tegration of the time-dependent Schrdinger equation for 
the level crossing models studied in the present paper 
(most of which are very slowly convergent from a numer- 
ical viewpoint) becomes increasingly difficult and time 
consuming as we approach this limit. In Appendix |a|, 
we present the numerical approach used in our studies, 
which has been specifically developed for highly accu- 
rate numerical integration of slowly convergent two-state 
problems, such as the present ones. The need for sofisti- 
catcd numerical approaches emphasizes the usefulness of 
the DDP formulas (^) and ([lO]), which become increas- 
ingly accurate as the adiabatic limit is approached. 



both a and 7 being positive. Then Eqs. ( p^ ) become 

0(r)=a, A(t)=t(1± 7 V + ...), (14) 

where Q(r) = Q(t)//3, A(r) = A(i)//3, and the plus (mi- 
nus) sign is for the superlinear (sublinear) case. 

The direct treatment of model (14) involves a third- 
order algebraic equation which is too cumbersome. More- 
over, for a negative nonlinearity, model ( |l4| ) involves two 
unwanted additional spurious crossings. We avoid these 
drawbacks of model (|14j) by using instead two other mod- 
els which contain positive and negative cubic nonlincar- 
ities, have a single level crossing, involve dealing with 
quadratic equations, and allow simpler derivations. 

A. Superlinear model 

The first model we consider is 



O(r) = a, A(t) = iVl + 27 2 



(15) 

For 7 = 0, model ( |l5[ ) reduces to the LZ model. Near the 
crossing (r = 0), we have A(r) == t (l + 7V 2 + ...), i.e., 
model (|l^) reduces to model ( |l4| ) with a positive nonlin- 
ear term. Moreover, as t — ► ±00, the detuning diverges 
quadratically, i.e., faster than the (linear) LZ detuning. 
Hence, we call model ( |l5|) the superlinear model. 

This model has two transition points in the upper com- 
plex half-plane, 



1 a 
J 



where £ = 2^/2a^. For 7 — > 0, we have r c 



(16) 



and 



7" — I 

point. 



and we recover the (single) LZ transition 



III. PERTURB ATIVE NONLINEARITY 

We wish to estimate the probability for nonadiabatic 
transitions P in the case when the linear LZ detuning (^) 
is perturbed by a small cubic nonlinearity in the vicinity 
of the crossing (at t = 0), 

f2(t) = const, A(t) = f3 2 t + Tt 3 + .... (12) 

When r > 0, the detuning passes through resonance in 
a superlinear manner (|A(t)| > |A£z(i)|), whereas when 
r < 0, the detuning passes through resonance in a sub- 
linear maimer (|A(i)| < |Alz(^)I)- 11 is m ore convenient 
to work with dimensionless quantities and we choose (3 
to define our frequency and time scales. We define the 
dimensionless coupling a, the nonlinearity coefficient 7, 
and the time r as 

a =7;, 7= Mr. r = pt, (13) 



1. The case of small aj 

For £ = 2V2a7 < 1, both transition points are purely 
imaginary, with Imr~ < Imr+. Hence, we can take the 
contribution from r~ only. We make the substitutions 
t = iy and r c 7 = iy c in Eq. (||) and obtain 

V (t~) =2i V" 2 - V 2 + 27V dy 
Jo 

= \nar-F ( i -A; 2; 1 ~ V* ~ 

2 V 2 l + v/W 2 / 

= i^(J,|;2;f 9 ) 1 (17) 

where F(a, 6; c; z) is the Gauss hypergeometric function 
@. The DDP formula gives 

P^ e — 2 ^(if;2;€ 2 ). ( 18 ) 



3 



For £ = 2y/2aj < 1, Eq. Q reduces to 



P 



3 - ff a 3 (l+iot a 7 a +-) 



(19) 



Comparison with the LZ formula (g) shows that the 
superlinearity reduces the probability for nonadiabatic 
transitions. 

The DDP result ( [l9|) is formally valid in the near- 
adiabatic regime (a» 1). Moreover, in the derivation of 
Eq. (19) we assumed that 2\/2aj <§; 1. Hence, we should 
have 



1 



> a > 1. 



(20) 



We have to account also for the fact that both the other 



zero point r+ and the singularity of model (15), situated 
at To = % /jV%, should be well separated from the transi- 
tion point t~ . It is easy to verify that in the present case 
(27207 < 1) we have Imr c < Imr+ < Imro and hence, 
we should only have Imr+ — IrriT" 3> 1. This is indeed 
the case if condition ( p(i| ) is satisfied. 

Condition ( p0| ) limits the applicability of Eq. (|lj]) to 
very small nonlinearities only. We have found, how- 
ever, through extensive comparison with numerical sim- 
ulations that this condition can be relaxed to 



1 



2V27 



> 



> 1. 



(21) 



2. The case of large 07 

For £ = 2v / 2a 7 > 1, the transition points ( |l6| ) are com- 
plex, with equal imaginary parts and real parts which are 
equal in magnitude but opposite in sign. Then, following 
Eq. ([HJ) , we have to take the contributions from both of 
them. The integrals T>(t^) are given by 



1.9. ±±wE 



1 



(22) 



where ±T> r and T>i are the real and imaginary parts. The 
factor s (pi] ) are T± = — 1 and the generalized DDP for- 
mula (10) gives 



P ~ 4e~ 2V * cos 2 V r 



(23) 



In Fig. the probability for nonadiabatic transitions 
for the supcrlinear model (|l5|) is plotted against the di- 
mensionless coupling a for nonlinearity 7 = 0.25. The 
DDP approximations ( |l9| ) and ( p3| ) are seen to fit the 
numerical results very well. As seen from the figure, this 
is the case even for small a, while the DDP formula is 
supposed to be valid for a > 1 only. This is a conse- 
quence of the fact that the DDP formula provides the 
exact result for the LZ model for any a. 
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FIG. 2. The probability for nonadiabatic transitions P 
for the superlinear model ( |l5| ) plotted against the dimension- 
less coupling a = Q//3 for 7 = 0.25. The numerical results 
are shown by dots while the DDP approximations ( p~9|) (for 
a ^ \/2) and (|23) (for a > y/2) by solid curves [the small-a 
approximation (jl9|) is almost indistinguishable from the nu- 
merical results] . The LZ probability ™ is shown by a dashed 



The condition of validity of Eq. (|23j) is expected to 
be £ = 2\/2Q!7 > 1. However, for large coupling a and 
nonlinearity 7, the singularity point to = «/7v / 2 gets 
closer to the real axis than the transition points which 
makes the DDP result ( p3| ) inaccurate for very large 07; 
indeed, we have verified this numerically (not shown in 
Fig. ||) . We should emphasize that the case of large non- 
linearity resembles more that of the essential nonlinearity 
(Sec. |yp and hence, it is not interesting in the context 
of perturbative nonlinearity, considered in this section. 



B. Sublinear model 



The sublinear model is defined as 



fi(r) = a, A(t) 



</l + 4 7 V 2 ' 



(24) 



For 7 = 0, it reduces to the LZ model (|4|). Near the cross- 
ing the detuning behaves as A(r) = r (l — 7 2 r 2 + ...), 
i.e., model (^4|) reduces to model ( [Til) with a negative 
nonlinear term. As r — > ±00, the detuning diverges as 
|r|2, i.e., slower than the (linear) LZ detuning. Hence, 
the name sublinear model. 

This model has a single transition point in the upper 
half plane given by 



t c = iaJ V4c? 7 4 + 1 - 2a 2 7 2 
The integral T> (r c ) is given by 



(25) 
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FIG. 3. The probability for nonadiabatic transitions P for 
the subhnear model (^) plotted against the dimensionless 
coupling a = fl/0 for 7 = 0.25. The numerical results are 
shown by dots and the analytic approximation (E^) by a solid 
curve. The LZ probability M is shown by a dashed curve. 
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FIG. 4. The probability for nonadiabatic transitions P 
plotted against the dimensionless nonlinearity coefficient 7 
for the superlinear model (jl^) (full circles) and the sublinear 
model ( j24[ ) (open squares). The coupling is a = ft/0 = 1. 
The solid curves depict the analytic approximation ( J30[ ) de- 
rived from the DDP formula. 



V ( Tc ) = ij, 



where 



J = 2 




V 1 - 4 7 2 r 



dy. 



(26) 



(27) 



where the substitutions r = iy and r c = iy c have been 
made. The DDP formula (g) gives 



-2.7 



(28) 



In the limits of small and large 07, we find 

e -™ 2 (i-f«V-H, (v^o7<l), (29a) 

(29b) 



P^ e - 2 «/7+V(16a 7 3 ) + - ; (V2q, 7 >1). 



Comparison with the LZ formula (||) shows that the sub- 
linearity increases the probability for nonadiabatic tran- 
sitions. Note also that due to the absence of other tran- 
sition points in the upper half-plane, there are no oscil- 
lations in P. 

In Fig. the probability for nonadiabatic transitions 
P for the sublinear model (HJ) is plotted against the di- 
mensionless coupling a for nonlinearity 7 = 0.25. The 
DDP approximation ( p8|) is seen to fit the numerical re- 
sults very well. 

The condition of validity of the DDP approximation 
( p8| ) is expected to be a > 1. In fact, as Fig. || shows, 
Eq. (|28| ) is quite accurate for a < 1 too, which is related 
to the fact the DDP formula provides the exact LZ prob- 
ability, as noted above. On the other hand, we have to 
account for the presence of a singularity at tq — i/2"/. 
The relation Imr c < Imrg is always fulfilled which means 



that the DDP formula ( |28| ) should be accurate if the two 
points t c and r are well separated. For large cry, how- 
ever, they approach each other [see Eq. (p5|)], which ex- 
plains the small inaccuracy of Eq. ( pl| ) seen in Fig. ^. 



C. Discussion 



Comparison of Eqs. ([l4j), ( fi9| ) and (29a) shows that 
the probability for nonadiabatic transitions for model 
(|l2|) for small nonlinearity is given by 



P 



P LZ e--^ T fP\ 



(30) 



This equation is valid for 2 r//3 6 < 1. Since the devia- 
tion from the LZ probability Plz depends exponentially 
on the coupling f2 and the nonlinearity coefficient T, this 
deviation can be very large. On the other hand, how- 
ever, this difference emerges when f2//3 > 1; then P is 
very small (virtually zero), which explains why the LZ 
formula is often found to be more accurate than antici- 
pated. 

In Fig. |I| the probability for nonadiabatic transitions 
P is plotted against the nonlinearity coefficient 7 for the 
superlinear model (|l^) and the sublinear model (p4j) . For 
7 = 0, both the superlinear and sublinear probabilities 
are equal to the LZ value Plz ~ 0.043. As 7 departs 
from zero, the superlinear probability decreases while the 
sublinear probability increases, in agreement with our 
analytic result (p0|). This qualitatively different behav- 
ior is readily explained by looking at the quasienergies 
±£(t). For the superlinear model (|l5|), the avoided cross- 
ing between the quasienergies is sharper than in the LZ 



5 



case, while for the sublinear model (p4[), it is flatter than 
the LZ one (see Fig. |l|). Consequently, the "duration" 
T of the avoided crossing for the superlinear model is 
shorter than for the LZ model, while for the sublinear 
model it is larger, T super < Tlz < %ub, which means 
that the same relation should apply for the probabilities, 



P. 



super 



< Plz < P, 



sub • 



D. Comparison with the Allen-Eberly-Hioe model 

An argument in favor of Eq. ( 30T) can be derived from 
the Allen-Eberly-Hioe model |34B3], 



Q(t) 



A 



cosh(i/T) ' 



A(t)=Btsnh(t/T), (31) 



which is a particular case of the Demkov-Kunikc model 
jMp2|. Model (31) can be solved analytically and gives 



P = 



cosh 2 (WB 2 - A 2 T) 
cosh 2 (7rBT) 



(32) 



It has been shown in |56| that this model is a member of 
a class of infinite number of models, for all of which the 
nonadiabatic probability is given by Eq. ([32]). Another 
member of this class is the model 



Q(t) = A, A(t) =Btsm(t/T), 
(-±ttT ^ t £ ±ttT). 



(33) 



We can make model ( |33|) behave like the perturbatively 
nonlinear model ( |l2] ) (i.e., make the coupling duration 
infinite and the detuning divergent) by letting T — > oo 
and B — > oo, while maintaing the slope at the crossing 
B/T constant. Then Eq. © reduces to 



exp 



A 2 



A 2 T 

' 4B 2 



(34) 



The same result follows from the DDP approximation 
( |30| ) by accounting for the fact that for model (|33|), a — 
A^/TjB, (3 2 = B/T, and T = B/3T 3 . 



IV. ESSENTIAL NONLINEARITY 

The essentially nonlinear model is defined by 

fi(t) = const, A(t) = (3 N+1 t N , (35) 

where N = 3,5,7,... is an odd number. Model (|35| ) 
cannot be linearized in any vicinity of the crossing. 

The zero points of the quasienergies @ in the upper 
t-plane are 

r k c = a i/" e *r(2fc-i)/(3A0 (fc = l, 2 ,...,A), (36) 



where r = (3t and a = f2//3. With the exception of 
Tc + 2 , which is imaginary, the transition points are 
grouped in pairs that have the same imaginary parts but 
opposite real parts. The significant difference between 
the LZ model and the essentially nonlinear model is that 
^_(N+i)/2 j g ^ e on jy transition point in the former, while 
in the latter it is the one farthest from the real axis and 
hence, with the smallest contribution to the nonadia- 
batic probability. The largest contribution comes from 
the points t\ and which are closest to the real axis. 



The integrals 2?(t*) are given by 



V{r k c ) = 2v N a 



(N+l)/N ott(2A-1)/(2JV) 

t- 1 



(37) 



where the number v-^ denotes the value of the integral 

(38) 



v N = [ \/\ - x 2N dx. 
Jo 



The first few values of are v\ — 7r/4 w 0.785, ^3 w 
0.911, ^5 w 0.944, v 7 w 0.959. The factors (|ll|) are given 
by Tfe = (— l) fc . The generalized DDP formula ( |Io| ) gives 



2 ( — l) c~ ~2N~ COS f Tj COS — ^~7T — 7f 



fc=l 

+ (_ 1 )(A r +D/2 e - t , 



2A 



(39) 



with 77 = 2z//va( JV+1 - ) / Ar . For large a, the main contri- 
bution to P comes from the first term (k = 1) in the 
sum. 

The probability (|3^) is plotted in Fig. || against the 
coupling a and compared with the exact numerical re- 
sults for N — 3 (top figure), N = 5 (middle figure), and 
N = 7 (bottom figure). Obviously, the probability ( |39| ) 
is qualitatively different from the LZ probability Plz- 
First, there are oscillations which appear due to the ex- 
istence of multiple transition points. From another point 
of view, the oscillations appear because the nonadiabatic 
coupling d(t) has two peaks, in contrast to the single- 
peaked nonadiabatic coupling in the LZ model. Second, 
the transition probability is much larger than Plz which 
can be explained by the fact that in the quasienergy pic- 
ture, the avoided crossing has a much longer duration 
(see Fig. |). 

These differences can be verified experimentally. For 
example, an essentially nonlinear crossing arises when a 
two-state atom is excited by a frequency modulated laser 
pulse with a supergaussian time dependence [i.e., Cl(t) oc 
e -(*/ T ) ] anc | the frequency modulation is produced by 
the self-phase modulation technique, in which the phase 
shift (pit) is proportional to either the amplitude of the 
field [i.e., (pit) oc f2(t)] or the intensity of the field [i.e., 
(pit) oc ^ 2 (i)] (37]]. Since the detuning is proportional to 
(pit), we shall have A(t) oc i 2 ™" 1 with a crossing at t = 0. 
As the hypergaussians are almost constant near t = 0, 
this example reproduces model (|35|) very well. 
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FIG. 5. The probability for nonadiabatic transitions P for 
the essentially nonlinear crossing model ( |35| ) plotted against 
the dimensionless coupling a = £1/(3 for iV = 3 (top figure), 
N = 5 (middle figure), and N = 7 (bottom figure). In each 
figure, the numerical results are shown by dots, the analytic 
approximation ( |39| ) by a solid curve, and the LZ probability 
(h) by a dashed curve. 



time dependence near the crossing. In the latter, the de- 
tuning is proportional to an odd power of time, A cx t N 
with N = 3, 5, 7, . . .. 

For the perturbative nonlinearities, the probability for 
nonadiabatic transitions P is larger than Plz for sub- 
linear nonlinearity and smaller than P^z for superlin- 
ear nonlinearity. For both the superlinear and sublinear 
cases, the appreciable deviations from Plz emerge for 
large coupling when P is already very small, virtually 
unobservable. This fact explains why the LZ model has 
often been found to provide more accurate results than 
anticipated. We have provided a simple analytic estimate 
for the deviation as a function of the nonlinearity, which 
should be a useful criterion for estimating the applicabil- 
ity of the LZ model to any level crossing case, as far as 
perturbative nonlinearity is concerned. 

For the essential nonlinearity, we have found that the 
nonadiabatic probability P is both quantitatively and 
qualitatively different from Plz because on the one hand, 
it vanishes in an oscillatory manner as the coupling in- 
creases, and on the other, it is much larger than Plz- 
From a mathematical viewpoint, these differences de- 
rive from the existence of two complex transition points 
equally close to the real axis, both being closer than the 
single LZ transition point. From a physical viewpoint, 
the differences can be explained by the fact that the 
avoided crossing is flatter and of longer duration than the 
LZ one and the nonadiabatic coupling has two peaks, in 
contrast to the single-peaked LZ nonadiabatic coupling. 
We have suggested an experimental situation when this 
deviation can be observed. 

We have limited our analysis to the case when the de- 
tuning is an odd function of time [Eq. (|^)], which excludes 
asymmetric level crossings, e.g., quadratic corrections in 
the perturbative models. It would be interesting to es- 
timate the effect of such quadratic terms although then 
the DDP treatment is more complicated. It should be 
noted that the parabolic level crossing model considered 
in f30|l involves either two level crossings or no crossing 
and hence, it is different from the present case of a single 
level crossing. 

Finally, we point out that the detuning nonlinearity is 
only one of the sources of possible inaccuracies in applica- 
tions of the LZ model. Others include the finite coupling 
duration Bql and the nonzero transition times [BSfl . 



V. SUMMARY AND CONCLUSIONS 

We have examined the effect of nonlinearity in the de- 
tuning at the level crossing on the probability for nona- 
diabatic transitions P. Our analysis has been based 
upon both analytic approximations derived by using the 
Dykhne-Davis-Pechukas formula and numerical calcula- 
tions. We have distinguished two types of nonlinearities: 
perturbative and essential. In the former type, the non- 
linearity appears as a correction to a dominant linear 
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APPENDIX A: NUMERICAL INTEGRATION OF 
SLOWLY CONVERGENT TWO-STATE 
PROBLEMS 

In this Appendix, we describe the numerical algorithm 
we have used to integrate the Schrdinger equation ([!]). 
The Landau-Zener model is notoriously known for its 
slow numerical convergence because the amplitude of the 
oscillations, which appear in the time evolution of the 
transition probability P, only vanishes as |f| _1 in the 
diabatic basis |3(|. This is a consequence of the slow 
(linear) divergence of the detuning. Of course, the final 
LZ transition probability (Q) is known exactly, but nu- 
merical integration is necessary when the time evolution 
is needed. The situation is even worse for the sublinear 
model ( p4| ) in which the oscillations amplitude vanishes 
only as |f|~ 5 . One of the possibilities to alleviate this 
problem is to carry out the numerical integration in the 
adiabatic basis, where the oscillation amplitude vanishes 



as If I for the LZ model 



and as If I for the sub- 



linear model. This improvement may be insufficient if 
the transition probability is very small, as in the present 
paper. Yet another problem arises from the finite initial 
(large negative) time, because it introduces additional 
oscillations in P |3(|. For the LZ model these problems 
can be resolved by starting the integration at the crossing 
(f = 0) and propagating the solution towards the desired 
(positive or negative) time |56|,[38| . The initial conditions 
in this approach require the values of P and a few of its 
derivatives at f = 0, which can be found exactly p6| , p8[ . 
Unfortunately, this approach does not apply to the mod- 
els in the present paper. 

We have used a numerical approach which combines 
and generalizes ideas by Bambini and Lindberg and 
Berry and Lim ]4C|] . The approach is based on two con- 
cepts. First, we perform the numerical integration in a 
superadiabatic (SA) basis ]4C|] , rather than in the usual 
diabatic or adiabatic bases. The successive SA bases are 
obtained iteratively J40|-[f2| and they are not the more fa- 
miliar superadiabatic bases obtained as truncated asymp- 
totic series in the adiabatic parameter 43 1. The n-th 
order SA states arc defined as the instantaneous normal- 
ized eigenstates of the Hamiltonian in the SA basis of 
order n — 1. For instance, by diagonalizing the Hamilto- 
nian in the diabatic basis (which is the SA basis of order 
n = 0) we obtain the adiabatic basis (SA basis of order 
n = 1). The recursive relations between the "couplings" 
f2„(t) and the "detunings" A„(f) in two successive SA 
bases are given by 



A„ +1 (f) = y/Ql(t) + Al(t), 

_ , . o„(t)A„(f) - n„(f)A„(f) 



(Ala) 
(Alb) 



Here A (f) = A(f) and Qo(t) = f2(i) are the detuning and 
the coupling in the diabatic basis, whereas Ai(f) = £ (f) 
and f2i(f) = i?(f) are the quasienergy and the nonadi- 
abatic coupling in the adiabatic basis. As far as the 



level crossing models in the present paper are concerned, 
the advantage of using the n-th SA basis is that the os- 
cillations in the transition probability evolution, whose 
amplitude is proportional to the ratio Q n (t) / A n (t) at 
large times, vanish much more quickly. This is so be- 
cause at large times, the SA "detuning" A„(f) diverges 
in the same manner as A(f), while the SA "coupling" 
fi n (f) vanishes as |fA(f)| _Tl . Hence, the oscillation am- 
plitude vanishes as |f|~™|A(t)|~ n_1 . We have used the 
third SA basis (n = 3), in which the oscillation ampli- 
tude vanishes as |f| -7 for the LZ model ([!]), as |f| -5 for 
the sublinear model (p4|), as |f|~ n for the super linear 
model ( |l5|) , and as Ifl^^^ -3 for the essentially nonlinear 
models (pq). Moreover, this approach provides the possi- 
bility to check the accuracy by calculating P in different 
SA bases. 

As we have pointed out above, another problem, which 
cannot be resolved merely by the choice of basis, is the 
finite initial time. We have overcome it in the manner of 
Bambini and Lindberg p9| by using the symmetry of the 
two-state problem in the n-th SA basis. The Bambini- 
Lindberg approach is based on a connection between the 
two-state evolution matrix U(f, 0), describing the evolu- 
tion from time f = to time f, and the evolution matrix 
U(f, — f), describing the evolution from — f to f. This 
approach allows to start the integration at f = 0, propa- 
gate it towards +00, and stop the integration when some 
convergence criterion is fulfilled [for instance, we have re- 
quired that the oscillation amplitude at time f is smaller 
than 10~ 4 P(f)]. It is easy to see that this approach is 
much faster than merely twice compared to the standard 
one (starting at large negative time), in which a conver- 
gence check would require to start the entire integration 
again at a larger negative time. The Bambini-Lindberg 
approach needs to be generalized because it applies to 
the case of a coupling and a detuning that are both even 
functions of time, while in the present case, the n-th SA 
"detuning" is an odd function in the diabatic basis and 
an even function in the SA bases, while the n-th SA "cou- 
pling" can be either an even function (for odd n) or an 
odd function (for even n). The relations between U(f, — f) 
and U(f, 0) for the four possible combinations of sym- 
metries in the detuning and the coupling (odd or even 
functions) have been derived in E3] in a similar manner 
as in |39|. We present the results in Table |. In the table, 
b\ = &i(+oo) and &2 = &2(+oo) are the so-called fun- 
damental solutions in the corresponding n-th SA basis, 
i.e., the solutions for the initial conditions fei(0) = 1 and 
&a(0) = at f = 0. As we can see from the table, in the 
different SA bases, the transition probability is expressed 
in terms of the fundamental solutions differently. More- 
over, depending on the SA basis, the desired probability 
for nonadiabatic transitions is equal to the probability 
of no transitions (for n = 0, 2, 4, . . .) or to the transition 
probability (for n = 1, 3, 5, . . .). 

The last case in Table |, in which both the "coupling" 
and the "detuning" are odd functions of time and which 
is given for the sake of completeness, is not interesting 
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in the context of the present paper but it represents an 
interesting effect called symmetry forbidden transitions, 
in which the system returns to its initial state in the end 
of the interaction |^a| . This return is only determined by 
the symmetry of the Hamiltonian and does not depend 
on its particular details. Moreover, it remains valid in 
the general iV-state case as well. 



Symmetry 


U„(+oo, -oo) 


Cases 


P 


n n {-t) = n n {t) 

An(-t) = A„(t) 


U X1 = b{ + bf 
U 12 = 2flm(& 1 & 2 ) 


odd n 


Pl2\ 2 


n„(-t) = Q n (t) 

A n (-t) = -A„(t) 


t/n = \bi\ z -\b 2 \ z 
U 12 = -2Wb* 2 


n = 




n n (-t) = -a n (t) 

A n (-t) = A n (t) 


Un =b{- bf 
U 12 = -2Re(6 1 6 2 ) 


even n 




n n (-t) = -n n (t) 

A n {-t) = -A n (t) 


Uu = l 

u 12 = o 







TABLE I. The relations between the fundamental solu- 
tions 61 (+00) and 62 (+00) in the n-th superadiabatic basis 
and the probability P for nonadiabatic transitions. In the el- 
ements of the evolution matrix U n (+oo, —00) in the n-th SA 
basis as well as in the fundamental solutions we have omitted 
for simplicity the arguments ±00 and the label n. 
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